Fractional quantum Hall physics in rotating ultracold Rydberg gases 
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We study ultracold Rydberg-dressed Bose gases subject to artificial gauge fields in the fractional 
quantum Hall (FQH) regime. The characteristics of the Rydberg interaction gives rise to interesting 
many-body ground states different from standard FQH physics in the lowest Landau level (LLL). 
The non-local but rapidly decreasing interaction potential favors crystalline ground states for very 
dilute systems. While a simple Wigner crystal becomes energetically favorable compared to the 
Laughlin liquid for filling fractions v < 1/12, a correlated crystal of composite particles emerges 
already for v < 1/6. The presence of a new length scale, the Rydberg blockade radius ob, gives rise 
to a bubble crystal phase when the average particle distance becomes less than as and v < 1/4. 
For larger fillings indications for strongly correlated cluster liquids are found. 



In the context of the fractional quantum Hall effect 
(FQHE) interesting many body ground states exists, 
some of which carry fractional and non-abelian braid- 
ing statistics [1-5 . This makes them ideal candidates for 
a topological quantum computer [6j [7] and gives access 
to topological quantum phases [8 . For an experimen- 
tal study a high degree of control is required and quasi- 
2D ultracold gases of atoms in artificial magnetic fields 
have been suggested as potential candidates [9 . Effective 
magnetic fields are generated either by rotation [lOj [11] 
or by using light-induced gauge potentials [T2UT4] . When 
a weakly interacting superfluid Bose gas is rotated suffi- 
ciently fast an ordered vortex-lattice forms p~5l [16] and 
the lowest Landau level (LLL) regime can be reached [T7] . 
Although there is no integer quantum Hall effect based 
on the Pauli principle for bosonic atoms, repulsive inter- 
actions can give rise to highly correlated FQH states for 
fillings v < 6 [9 when only the LLL is occupied [T8H20] . 
However, despite the experimental progress strongly cor- 
related phases have not been realized yet, which can par- 
tially be attributed to the rather small interaction ener- 
gies in atomic gases. In the present paper we propose an 
alternative using 1/r 6 van-der-Waals (vdW) interactions 
in Rydberg states. The associated energies can be orders 
of magnitude larger than those achievable with contact 
or magnetic dipole-dipole interactions. We theoretically 
analyze the phase diagram in the LLL using exact diag- 
onalization (ED) as well as variational methods. Similar 
to dipolar gases with 1/r 3 interaction, which have been 
extensively studied in the past [21-23 , we find a transi- 
tion from Laughlin (LN) liquids [24] to crystalline ground 
states for very dilute systems, however with an extended 
filling region where only composite particles can crystal- 
lize at zero temperature. Additionally the presence of 
a new length scale, the Rydberg blockade radius gives 
rise to a clustering mechanism supporting bubble crystal 
ground states. 

Recently there has been a lot of interest in atoms ex- 
cited to high Rydberg states and there has been con- 
siderable experimental progress to make these systems 
accessible [25H3T] . Of particular interest in the present 
context are atoms excited by far-off resonant laser radia- 



tion. In this case, called Rydberg dressing (see figjT](a)), 
the atoms essentially remain in their ground state but 
show an effective interaction 1321 l33l 
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For large particle separations r the interaction poten- 
tial is of vdW type - r~ 6 , where C 6 = (n/2A) 4 C 6 
describes the effective interaction strength of Rydberg 
dressed atoms with Cq being the bare vdW coefficient, 
Q <C |A| the Rabi- frequency of the laser excitation and 
A its detuning. Since bare Rydberg interaction energies 
are typically many orders of magnitude stronger than 
e.g. magnetic dipole-dipole interactions, the excitation 
gap of FQH states can easily reach the energy separation 
between Landau levels (LLs, see supplementary mate- 
rial for specific data). Most importantly the interaction 
potential ([!]) flattens off below the blockade radius a#, 
thereby defining a new characteristic length scale. 




FIG. 1. (color online), (a) An ultra-cold Bose gas is dressed 
with Rydberg-excitations \r) by illuminating with a far off- 
resonant (detuning A) laser beam Q. Rotating the gas with 
frequency cjo and applying a weak harmonic trap u± in ra- 
dial and a tight one in the transverse direction (uj z ^> u±) 
the quasi-2D FQH regime can be reached, (b) Contour plot 
of pseudopotentials Vm for the potential ([T]), where m was 
treated as continuous parameter for better illustration. Black 
solid lines denote contours of constant (continuous) pseudopo- 
tential. 
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When a Rydberg-dressed Bose gas is set into rota- 
tion with angular frequency uo in a radial trap of fre- 
quency uj±, the corresponding length scale is given by 
£ c = (2m a ^/ujQ + uJ^jK) -1 ! 2 [34 with m a the atomic 
mass. As will be shown the competition of the two length 
scales £ c and as leads to interesting new physics. Already 
in the mean-field regime where the rotation frequency is 
well below the deconfinement limit (uj± > ujo) there is 
e.g. a competition between a vortex-lattice and a super- 
solid phase [35]. Here we are interested in the regime of 
strong correlations. To this end we consider a Rydberg- 
dressed Bose gas in the LLL regime for fillings v < 1, 
assume a quasi-2D gas, neglecting finite-thickness effects 
and use uo±_ = uoq (£ c —> £b- magnetic length). For sim- 
plicity we consider only even values of \ jv corresponding 
to bosonic LN states. 

The interaction projected to the LLL is described by 
Haldane's two-particle pseudopotentials V m [36] . where 
the integer m denotes the relative angular momentum. 
For the interaction ([!]) they are plotted in fig[]Jb). At 
given as the first V m are approximately equal until the 
radial extend R m = \j2m£B of the wavefunction in the 
relative coordinate is ~ a#. The subsequent pseudopo- 
tentials decrease quickly. Specific numbers are given in 
the supplementary material. The two regions, separated 
by the curve \Jlm£B = clb give rise to qualitatively dif- 
ferent physics. We first consider the case of a small block- 
ade radius, where its effect can be disregarded. 

The v = 1/2 and v = 1/4 ground states. For 
o^b/£b the first two pseudopotentials diverge 

y °"8^UJ ' ^(°- 571 - ln ^)^ (2) 

while all V m> 2 converge. Therefore, the v = 1/2,1/4 
bosonic LN states are exact ground states. ED calcula- 
tions in spherical geometry for small particle numbers TV 
confirm this for small but finite values of a^/^: even 
for clb/£b — 2 we find large overlaps squared to the LN 
states, being 0.992 for v = 1/2 (N = 10) and 0.974 for 
v = 1/4 (N = 7). 

Ground states at small fillings v < 1/4. Now we want 
to address the physics for small fillings, i.e. v < 1/4, 
again for a# < £b- Due to the non-local interaction 
a natural ground state candidate in this regime is the 
noncorrelated Wigner crystal (NWC) [37] described by 
the wavefunction 

VnwcW = S jj e -(ki-^l a +* J fl;-*;H J )/4 > (3) 
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where Rj G C defines a lattice in the plane of complex 
coordinates Zj = (xj + iyj)/£B and S stands for complete 
symmetrization. To see if the NWC could be a ground 
state we compare the corresponding variational energy 
per particle £nwc to that of the LN state. To this end 
we generalize the expression for 6nwc derived in [37] to 
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FIG. 2. (color online). Comparison of the variational energies 
per particle eo (at fillings v — 1/2, 1/4, 1/6, ... corresponding 
to the bosonic Laughlin states) for the different trial wave- 
functions at ob = 0.2£b (LN - Laughlin states, NWC - non- 
correlated Wigner crystal of bare particles.) 



bosons: 

1(7 6 e -\Rj\ 2 /i 

6NWC = ^ ^ 1+ l e -\ Rj \*/2 ' K ^ = l^'D' ( 4 ) 

where Io,Jq denote Bessel functions. The integral can 
be handled numerically and the lattice sum can easily 
be performed. Like for electrons [37j [38, and dipolar 
fermions [39], the hexagonal lattice minimizes the NWC 
energy, see fig [2] 

In fig|2] the energies of different NWCs and the LN 
liquid are plotted for small as = 0.2£b and even values of 
\jv. The latter was calculated using the standard plasma 
analogy [24] for N = 100 particles for all examined v with 
Metropolis Monte Carlo (MC) [40 . We find a transition 
from LN to a NWC at v = ^nwc = 1/14, although for 
v = 1/12 the energy difference is so small that we can 
not exclude a transition here. This critical filling is far 
below the value of v « 1/7 for pure Coulomb and dipolar 
interacting fermions [38, 39| [4Tj . 

On the other hand a simple stability analysis shows 
that the lattice is expected to melt only for fillings 
v > ^Ld = 1/4. Specifically we calculated the Linde- 
mann parameter 7 = (Su 2 )/a from the phonon spec- 
trum in harmonic and nearest-neighbor approximation 
following [42 , where a denotes the (hexagonal) lattice 
constant and Su the atomic displacement from the lat- 
tice {Rj}. The phenomenological Lindemann criterion 
states that the crystal melts when 7 exceeds a critical 
value 7 C . In our case 7 = 0.57^ for as < £b at zero 
temperature, and j c « 0.28 [39j [43]. Thus the crys- 
tal is expected to melt only for v > 1/4, similar to the 
result found in dipolar systems [22 , 39 . This value dif- 
fers from the transition point to a NWC, z^nwc = 1/14, 
found previously. Instead of a NWC the ground state 
in the region z^d < v < ^nwc can however be a crys- 
tal of composite particles (CPs) [44], termed correlated 
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FIG. 3. (color online), (a) Variational energy per particle of 
the hexagonal CWC of^CP, ecwcW, in units of the LN state 
energy cln(^) at filling v. Results were obtained in Metropolis 
MC simulations with N = 91 particles, (b) Second order 
correlation function of the ground state for vdW interactions 
for N = 8 solid (N = 6 dashed) at v = 1/2 [y = 1/4). For 
large as symmetry is broken (L ^ 0, thick lines), in contrast 
to small as (L = 0, thin lines). 



FIG. 4. (color online). Phase diagram (for 1/v £ 2N) ob- 
tained from the comparison of variational energies of LN state 
(N = 100 using Metropolis MC) and bubble crystal. The op- 
timal number A; opt of particles per cluster (color code) was 
used in the BC region. The squares with error bars show 
the phase transition from LN to a symmetry-breaking state 
obtained numerically (ED, spherical geometry). 



Wigner crystal (CWC) [45j|46 . As for the LN states the 
first pseudopotentials are screened for CPs and only the 
effects of the long-range tails of the repulsive interaction 
potential remain. As the latter should be the same for 
CPs as for bare particles the stability arguments given 
above remain valid and a CWC should form already for 
v < ^Ld = 1/4. 

CWC variational wavefunctions of the form [6] 

^CWcW = ^LLL Y[(Zj ~ ^)^NWC(^) 
j<k 

were investigated in [45] using plasma analogy based MC 
simulations, and the critical filling v « 1/7 in electronic 
systems was correctly predicted. Here 7\ll is a projector 
to the LLL and (zj — z^)^ are Jastrow factors describing 
the formation of fi- composite bosons M CB (fermions^CF) 
for fi even (odd). We here performed similar MC simula- 
tions and calculated the variational energy of hexagonal 
CWCs for different \i. We take into account only direct 
but no exchange terms and check that the first exchange 
corrections are negligible. Our results are shown in figj3] 
(a) and we observe that for v < 1/6 CWCs are energet- 
ically favorable compared to LN states. We thus expect 
crystalline order already for v < z^d = 1/4 in agree- 
ment with the stability analysis and long before the sim- 
ple NWC becomes energetically favorable at v = z^nwc- 
More specifically we expect a 4 CB CWC ground state at 
v = 1/6 and v = 1/8. For v = 1/10 the 4 CB and 5 CF 
CWCs are equal in energy within our MC errors and we 
can not make a final conclusion about the underlying 
CPs. 

Our conclusion of crystallization for v < 1/6 is further 
supported by numerical results for small systems. Using 
ED in spherical geometry we obtained LN-like states at 
aB = 0.2£b for all v at the accessible particle numbers TV, 
however extrapolation of the gap to collective excitations 
yields a negative TDL value for v < 1/6. Additionally 
the second order correlations show strongly enhanced os- 



cillations as compared to those expected for LN liquids, 
while for small r the correlations still decay according to 
a power law, e.g. g^ 2 \r) ~ r 11 for N = 6 at v = 1/8. 
This is a direct indication that the ground state physics of 
these systems is dominated by M CPs with /i between 4 and 
6, as predicted by our MC results presented above. More 
details about these numerical signatures can be found in 
the supplementary material. 

Bubble crystal at large blockade radii. In the follow- 
ing we discuss the effects arising from the competition 
of the magnetic length £b with the new length, the 
blockade radius a#. Specifically we consider the case 
where there are more than two particles per blockade 
area Ab — 7ra B ^ i.e. the regime > ^J~2JvIb- ED re- 
sults obtain ed on spheres suitable for LN states show that 
for as < \f2jvtB, the ground state is a LN state with 
vanishing angular momentum L = 0. For ~ y^/vls 
a transition to a L / state is observed, which after 
mapping from sphere (radius R) to plane (momentum 
fc), kR = L [48 , corresponds to a breaking of transla- 
tional invariance. Figj3] (b) shows numerical results for 
the two-particle correlation for v — 1/2 and v — 1/4 
for different ratios of as/^B- For small as (L = phase) 
we find LN-like correlations while for large as (L ^ 
phase) we find particle bunching at r = 0, which indi- 
cates clustering. Clustering of k particles can also be 
seen in the k + lst-order density correlations, 

g(k+l) {z) = ((^(O))**^)^) (*( )) fe ). 

E.g. for as = 2.9£#, k = 2 particles cluster, resulting in 
a very small # (3) (0) w 2 x 10~ 5 (^/2tt) 3 (for v = 1/2, N = 
8), while # (2) (0) = 0.9 (^/2tt) 2 is still large. 

These numerical findings can be explained as follows: 
For as > y/^/v^B where there are more than two par- 
ticles per Ab, the long-range contribution to the inter- 
action energy can be reduced by bringing the particles 
inside Ab together. At the same time there is no energy 
penalty from interactions inside Ab since the potential 
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([!]) is constant for r < as- This provides a pairing mech- 
anism. Assuming that all particles within Ab undergo 
clustering, one expects a transition from k — 1 to k par- 
ticles per cluster at d$ = y^2k/u • £b- 

By analogy we expect the formation of a NWC consist- 
ing of clusters, when the reduced cluster filling v c \ — vjk 
is small. This state is referred to as bubble crystal (BC) 
and was e.g. considered for electrons in a weak magnetic 
field where v > 1, i.e. beyond the LLL [49 . The present 
situation is comparable since the effective interaction in 
a higher LL is "smeared out" around r = 0. A BC phase 
was also predicted for dipolar bosons at filling v = 1/2 
with large finite-thickness effects [50]. To verify that BCs 
are indeed good ground state candidates for large block- 
ade radii we calculate the variational energy per particle 
by generalizing Q. This yields 

^bcO; v ) = k • e N wc = ^cl) + (fc - 1) Vb, (5) 

where the second term describes the binding energy per 
particle required for the cluster formation. In figj4] the 
ground state phase diagram is shown determined by com- 
paring the variational energies of the LN liquid, the NWC 
and the BC. Also shown are the transition points between 
symmetry conserving (L = 0) and symmetry breaking 
(L ^ 0) states obtained from ED. One recognizes that 
BC ground states with k > 2 exist for v < 1/4 and 
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FIG. 5. (color online). Qualitative form of the LLL phase 
diagram. 



Large filling and large as-' For v — 1/2 the BC is 
higher in energy than the LN state for all as and thus 
the nature of the ground state for large a# and large 
filling is an open question. Besides clustered crystalline 
states many different correlated cluster liquids have been 
proposed and discussed in the context of the FQHE in 
higher LLs. Our preliminary studies indicate that in this 
regime cluster liquid states with exotic properties emerge 
such as the Haffnian state (Hf) [51]. ED calculations for 
small systems (TV < 10) suggest that in our system the Hf 
is an incompressible ground state for a#/^# ~ 2.8. ..3.5, 
which is surprising as it is believed that it describes a 
compressible state [51]. The arguments used in [51] rely 
however on the fact that the effective composite interac- 
tion is repulsive for all r, and there are good arguments 
that this may not be the case for the present Rydberg 
system. Thus this parameter regime, which may be eas- 
ier accessible experimentally than the low-filling regimes 
discussed above, is potentially a very interesting one. In 
order to make more definite statements further studies 
are needed however. 

Summarizing, we have shown that Rydberg-dressed 
Bose gases can give rise to extremely strong interactions 
and a variety of interesting correlated phases not found in 
the standard FQH physics of the LLL. This has two rea- 
sons: the rapid fall-off of the interaction potential with 
distance and the competition of two length scales, a# and 
£b- A qualitative picture of the LLL phase diagram is 
shown in figj5j In the limit of pure van-der-Waals inter- 
actions (clb /@b 0) the v = 1/2, 1/4 Laughlin states are 
exact ground states. Although we mainly discussed filling 
fractions with even values of \ jv the entire region is de- 
noted by "FQHE" since we expect the standard bosonic 
FQH physics [20 to hold for all fractional fillings larger 
than v = 1/6. For v < 1/6 we find correlated Wigner 
crystal ground states, while a noncorrelated Wigner crys- 
tal ansatz predicts crystallization only at v = 1/12. For 
a B ^ \J~QvtB and v < 1/4 a transition to a bubble crys- 
tal is expected. Finally, for larger fillings and blockade 
radii the nature of the ground state is an open question, 
and for filling v = 1/2 we speculated on a connection to 
cluster liquid states in particular the Haffnian [51] . 

We thank J. Otterbach for helpful discussions and M. 
Baranov for drawing our attention to the crystal stability 
analysis. F.G. acknowledges financial support from the 
the Graduate School of Excellence MAINZ. 
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Supplementary Material for "Fractional quantum Hall physics in rotating ultracold 

Rydberg gases" 

F. Grusdt and M. Fleischhauer 



I. PSEUDOPOTENTIALS AND ENERGY SCALES 



In disc geometry the Haldane pseudopotentials are given by 

— / 

2 2m+l m! J Q 



V m = 



dr r 



l+2m 



V{r)e 



2 /4 



(6) 



where for bosons (fermions) only m = 0, 2, 4, ... (m = 1, 3, 5, ...) are relevant. In fig. 1(b) of the main text we plotted 
Vb,2,4,..( a #) together with contour lines of the continuous function V(m) = V m for m G M defined as in ([6| above. 
Table [T| list s some explicit values of V m . The characteristic energy scales for Rydberg- Rydberg interactions are limited 
by demanding LLL-approximation (i.e. Vq <C ftco> c ), different from all previously discussed interactions in cold gases. 



This is due to the large values Vq may take, as can been seen from tables [D] and [TTT] for Rb. There we assume a 
fixed magnetic length £b ~ l/^ m (corresponding to u c « 2tt • 130Hz which is a realistic value [1 j) as well as a fixed 
ratio ft/ A = 0.1. For the 6051/2 Rydberg state in Rb, C 6 (n = 60)/2tt = 0.14THz/im 6 and we used the scaling law 
Cq(ti) oc n 11 [2]. In table [n] A = 2ir • 1GHz is constant and the principle quantum number of the Rydberg state n is 
varied. In table III on the other hand, n = 46 is fixed and the detuning A is varied. Note that in the latter case it is 
even possible to tune the interaction potential adiabatically. 
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TABLE I. Leading order bosonic (m even) pseudopotentials V m for the potential V(r) = Cq/(t q + a%) from the main text. 
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TABLE II. Realistic blockade radii as and energies Vq for Rb with uj c w 2ty • 130Hz (i.e. £3 ~ 1/xm), A = 10Q = 2tv • 1GHz as 
a function of the principle quantum number n. 
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TABLE III. Realistic blockade radii as and energies Vq for Rb with uj c w 2ty • 130Hz (i.e. £b ~ 1/xm), n = 46 as a function of 
the detuning of the dressing laser A. Note that the Rabi frequency Q was chosen as Q, — 0.1 A. 



II. EXACT DIAGONALIZATION 

In addition to variational calculations we performed 
state-of-the-art exact diagonalization (ED) studies. To 



minimize the role of finite-size effects we did all our sys- 
tematic investigations in spherical geometry. We work on 
standard desktop computers and have excess up to the 
following particle numbers: 
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The dimension of the full Hilbertspace containing all an- 
gular momentum multiplets is given, since we do not ex- 
plicitly exploit that [H,L% ot ] = in the numerics. In the 
numerics the filling fraction v is adjusted by changing the 
magnetic flux N<$> through the surface of the sphere. For 
v = 1/2, 1/4, ... Laughlin states 



TV- 1 



We eliminate leading-order finite size effects by rescaling 
the magnetic length [3 , 



1% 



In simulations cle must be chosen small enough since 
when the blockade area fills half a sphere, i.e. when > 
7t^N<$>/8£b, we expect the formation of two clusters on 
opposite poles. We indeed found large overlaps to two- 
cluster trial wavefunctions in this case and no conclusions 
about the thermodynamic limit can be drawn. 

In disc geometry we only performed systematic studies 
for small cle < £b and found similar results as in spherical 
geometry. The accessible system sizes are the same as 
above, although requiring slightly more CPU time. 



III. NUMERICAL RESULTS AT v — 1/6, 1/8 

In the main text we summarized the most important 
numerical indications for crystallization at v = 1/6,1/8 
and for small < Here we give further details 

about these ED calculations in spherical geometry. We 
obtain substantial overlaps to Laughlin (LN) states for 
small systems, see table |IV| Besides a region of almost 
zero overlap (large a#) we find that for v = 1/6 LN 
states describe the ground state much worse than for v = 
1/2,1/4. In the main text we stated that for v = 1/6, 1/8 
finite-size extrapolations yield negative gaps to collective 
excitations. This can be seen in fig. [6] (a). The results for 
v = 1/4 are shown for comparison, and they clearly yield 
a positive gap in thermodynamic limit. The gap AE was 
obtained as the difference from ground to first excited 
state. The latter can well be described as a density wave 
and exhibits a roton minimum [4j [5] . Vanishing AE — » 
thus signals an instability with respect to crystallization. 

In fig. [6] (b) second order correlations are shown for 
v = 1/8 and = 0.2£#. As mentioned in the main 
text, we find strongly enhanced oscillations for N = 6 
signalizing a trend towards long-range order. For N = 5 
this effect is absent, which can be attributed to a more 
general finite-size effect: Already for the gaps (fig. [6] 
(a)) we observe a different behavior for N even/odd at 
v = 1/6,1/8. For N = 2,4,6, AE is smaller than for 



N = 3, 5 and we find the same behavior for the ground 
state energy. This is another indication for crystalliza- 
tion, since also for a classical vdW crystal on a sphere the 
ground state energy per particle shows these oscillations 
[? ]. The system is thus expected to be more suscep- 
tible to density modulations for even N. This in turn 
explains (i) the smaller gap and (ii) the enhanced oscilla- 
tions of the ground state correlations at N = 2, 4, 6. For 
v — 1/4 this finite-size effect is not observed, fully consis- 
tent with our result (see main text) that crystallization 
takes place only for v < 1/4. Finally we note that the 
discussed finite-size effect is due to incommensurability 
but not due to pairing since the latter would imply large 
g( 2 )(0) but small g^ 3 \0) which we do not observe. 
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FIG. 6. (color online), (a) Gap to collective excitations Ai? 
at v = 1/4, 1/6, 1/8 and for a B = 0.2£ B . Note that for v = 
1/4 energies were scaled down by a factor of 100. ED in 
spherical geometry was used and finite-size corrections were 
performed as described in sec. [Ii] Solid lines: quadratic fits, 
(b) Second order correlation functions of the ground state 
for van-der-Waals (vdW) interactions at as — 0.2£b and v — 
1/8, compared to those of the LN liquids at the corresponding 
system size. ED in spherical geometry was used. 
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TABLE IV. Overlaps squared of the ground state to the Laughlin states at fillings v — 1/2, 1/4, 1/6. Both were obtained using 
ED at N — 6 in spherical geometry. 
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